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Abstract 

The Spitzer Space Telescope provides a unique view of the Universe at infrared wave- 
lengths. Improved sensitivity and angular resolution over previous missions enable de- 
tailed studies of astrophysical objects, both in imaging and spectroscopic modes. Spitzer 
observations of active galactic nuclei can help shed light on the physical conditions of the 
central regions of these active glalaxies. 

The nearby radio galaxy Cygnus A is one of the most luminous radio sources in the 
local Universe. In addition to the high radio power, it is also very luminous in the infrared. 
New Spitzer spectroscopy and photometry of Cygnus A is combined with data from the 
literature at radio and sub-mm wavelengths. The resulting complication is modeled with 
a combination of: a synchrotron emitting jet, a burst of star formation, and emission from 
an AGN torus. 

The infrared emission in Cyngus A shows contributions from all three processes and 
the models are able to reproduce the observed emission over almost 5 dex in frequency. 
The bolometric AGN luminosity is found to be ~ 10 45 erg s _1 , with a clumpy torus size of 
~ 7 pc. Evidence is seen for a break in the synchrotron spectrum in the mid-infrared. The 
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relevant component of the infrared emission suggests Cygnus A has a star formation rate 
of ~ 20 M Q yr _1 . Even in the absence of the AGN, it would still be a luminous infrared 
source. 
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Chapter 1 

Introduction 



The advent of astronomical observations in non-optical wavelength regimes led to the 
discovery of bright celestial sources at radio wavelengths. Many of these sources of non- 
thermal radio emission have since been hypothesized to be powered by the release of 
gravitational energy from an accreting supermassive black hole at or near the center of a 
galaxy. These "active galactic nuclei" (AGN) represent some of the most energetic phe- 



nomenon in the Universe, and may play an important ro 



of galaxies and 



Nesvadba et al 



galax y clusters (e.g. 



Sanders et al. 



1988 



e in the formation and evo 



2008T ) 



Antonuccio-Delogu Silk 



ution 



2008 



There are many observational types of active galaxies, with drastically different mor- 
phologies and powers. Radio-loud AGN are classified as those sources whose radio flux 
exceeds the optical flux by at least a factor of ~ 10. The first radio-loud AGN to be 
discovered were associated with point sources at optical wavelengths. These "quasi-stellar 
objects" (quasars) are some of the most intrinsically luminous objects in the Universe. 

Frequently these radio- loud AGN exhibit extended structure in the radio wavelengths, 



and are classified as radio ga 



elliptical galaxies ([McCarthy 



axies. The host galaxies for these AGN are generally large 



1223). 
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There are multiple observational classifications of radio-loud AGN. Geometrical frame- 
works have be en suggested to unify s ome of these observational classes into a single type 



of object (e.g. lUrrv &; Padovani 



19951 ). Invoking both dust obscuration and source orien- 
tation, the appearance of a radio-loud AGN can be explained by the line of sight to the 
observer. 

Dust and gas obscuration predominantly affects wavelengths from the near infrared 
and shorter. Absorption from dust can obscure emission, especially from deeply embedded 
sources. At longer wavelengths the emission passes through the dust mostly unobscured, 
allowing direct studies of "obscured" regions. Additionally, the shorter wavelength radi- 
ation absorbed by the dust is re-radiated in the mid- and far-infrared. This makes the 
infrared the ideal regime to study bolometric output of an AGN. 

The radio galaxy Cygnus A is one of the brightest extragalactic radio sources visible 
in the sky. Due to its proximity, it can be studied at high spatial resolution, enabling 
a detailed study of the AGN and its environment. Studies of Cygnus A have indicated 
there is a powerful AGN at t he center, likely compa rable to quasars seen in the distant 



(therefore younger) Universe (jCarilli &: Barthel 



1996). Studying this nearby example can 



yield significant insight into the processes at work in distant quasars. 

Infrared measurements of Cygnus A can test the unification schemes by looking for 
the signature of an obscured AGN and comparing the observed infrared emission with 
models. Additionally, the total infrared emission from the AGN can be used to determine 
the bolometric luminosity, testing the idea that Cygnus A is a nearby example of a powerful 
quasar. 

Cygnus A also shows evidence of significant star formation around the nucleus. Dust 
enshrouded star formation has its primary signature in the infrared, lending further use 
to observations in this wavelength regime. 

The Spitzer Space Telescope (formerly the Space Infrared Telescope Facility) was 
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launched in 2003 and contains a unique suite of instruments for the mid- and far-infrared. 
With improved angular resolution and sensitivity over previous infrared satellites, it can 
provide significant additional detail over previous studies. Spitzer covers a wavelengt h 



200J). 



range of ~ 3 — 160 fim in a combination of imaging and spectroscopy (j Werner et aJj 

This work presents new mid-infrared spectroscopy of the nuclear regions of Cygnus A. 
Due to the significant energy output in this wavelength regime, interpreting the infrared 
emission from Cygnus A is critical to understanding the nature of the AGN and other 
nuclear activity. Three detailed, physically motivated models were combined to reproduce 
the observed spectral energy distribution covering 5 dex in frequency. 

The new data, plus data compiled from the literature, were modeled to explore the 
relative contributions from AGN activity and the circumnuclear star formation. The AGN 
activity is manifested as point-source emission and a relativistic jet. The jet dominates 
in the radio regime, while the nuclear point source emission is primarily evident in the 
infrared. The star formation is dust enshrouded so the resulting emission occurs primarily 
in the thermal infrared. 

Chapter [2] explores the motivation for this study by describing the observational zoo of 
AGN and the possible simplification which arises from a unification scheme. This chapter 
also describes the current body of knowledge regarding Cygnus A, which in turn provides 
the justification for the choice of models. The Spitzer Space Telescope and the instruments 
used to carry out new observations of Cygnus A are discussed in chapter [3l The modeling 
method and individual models used are described in chapter [U Finally, the results of this 
modeling and interpretation of the results are the subject of chapter 
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Chapter 2 

Active Galactic Nuclei and 
Cygnus A 



Active galactic nuclei (AGN) represent some of the most energetic phenomena in the 
Universe. They can influence the formation and growth of galaxies and galaxy clusters 
through feedback in the form of winds and jets. Likely powered by the extraction of 
gravitational energy from objects in the potential of a supermassive black hole (M. > 



1O°M , M Q = 1.99 x 10 33 g 



, they provide a laboratory to study high energy physical 



processes otherwise inaccessible (jBlandford fc Znajek 



1973). 



2.1 Active Galactic Nuclei 

The first (radio-loud) AGN were discovered via radio surveys in the 1950s. The relatively 

compact radio sources which were detected had no obvious optical counterparts. Follow 

up observations identified optical counterparts, ranging from optical point sources to giant 

elliptical galaxies. The optical spectra of these objects implied a large redshift, resulting 
lr rhe CGS system of units is utilized throughout this work, except where otherwise noted. 
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from a high recession velocity. The large recession velocities indicated these sources were 
very distant and therefore have high intrinsic luminosities. With efficienciesy of ~ 0.7%, 
nuclear fusion cannot provide sufficient energy output to power these objects. These 
objects are now believed to be powered by accretion of matter onto supermassive black 
holes (SMBH), where efficiencies approach 10%. 

There are a variety of observational types of AGN, broadly divided into two categories, 
"radio-loud" and "radio-quiet" , based on the observed ratio of radio to optical emission. 
Roughly 10% of AGN are "radio-loud" although this fraction varies with intrinsic power 
of the AGN. 

As the name implies, "radio-quiet" AGN are those with minimal radio emission, al- 
though they may show weak jets and radio lobes. Typically labeled "Seyfert Galaxies" 
(after Carl Seyfert), their optical spectra often show strong emission lines. In some cases, 
the emission lines are quite broad (~ few 1000s of km s _1 , "Seyfert I", "Type I" AGN). 
In sources with only narrow emission lines, they are still broad compared to a quiescent 
galaxy (~ few 100s of km a -1 , "Seyfert II", "Type II" AGN). 

"Radio- loud" AGN show strong radio emission. The radio morphology can vary dras- 
tically from object to object, from point source through large extended structures up to 
Mpc scales. Polarization studies of the radio emission, coupled with multiwavelength mea- 
surements of the flux, indicate the emission mechanism is synchrotron radiation from a 
population of relativistic charged particles moving in a magnetic field. "Radio-loud" AGN 
exhibit a similar range of optical spectra as their "radio-quiet" counterparts. As noted 

above, Cygnus A is a radio-loud AGN and its properties will be discussed in §2.21 

2 Efficiencies are defined as the fraction of the rest mass energy (E — mc 2 ) released as electromagnetic 
radiation. 
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2.1.1 Unified AGN Model 



There have 



aeen attempts to unify the zoo of observational AGN types into a consistent 



framework ( Urrv &i Padovani 



19951 ). Here the focus will be on unification of radio- loud 
AGN. The two major mechanisms used to explain the variety in observed AGN types are 
orientation-dependent obscuration and relativistic beaming of jets. 

Radio-loud AGN typically show extended structure at radio wavele ngths. The mor- 



pholo gy is generally divided into Fanaroff & Riley (FR) I and II objects (jFanaroff &; Rilev 



1974 ). FR I sources are brightest at the center, with emission tapering off away from the 
radio core. FR II sources are "edge brightened", showing distinct lobes with "hotspots" 
of radio emission towards the edges. It is widely believed that FR II radio sources contain 
relativistic jets, however the jets are often too faint to be seen. Figure 12.11 shows radio 
images of a typical FR I and FR II. 

The jet morphology of the FR II sources is evident, in figure ETT1 The overall source 
morphology is generally attributed to relativistic jets which propagate from the radio 
core with bulk relativistic motion. These jets likely terminate in the "hotspots" of radio 
emission as the jet runs into the tenuous intergalactic medium. Shocks at the interface 
likely re-accelerate the electrons and the radio plasma backflows to create the observed 
lobes of emission. 

As noted above there is considerable variation in optical spectra, with some objects 
showing narrow emission lines, and others with both broad and narrow emission lines. The 
"narrow" emission lines are still broader than lines typically seen in quiescent galaxies. 
Objects with broad emission lines are classified as "Type 1" AGN and those with only 
narrow lines as "Type 2". The variation in optical spectra is attributed to orientation- 
dependent obscuration of the AGN. 

Figure 12.21 shows a schematic diagram of the geometric structure invoked to unify 
radio-loud AGN. The central object (generally believed to be a supermassive black hole) 
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Figure 2.1: Left: 20 cm radio map of FR I radio galaxy 3c31 (jLaing et al.ll2008l ). Notice 
the emission tapers off away from the center of the source. Right: F R II source Cygnus 
A, with the hotspots and lobes clearly visible ( Carilli Barthel 19961 ). 



is surrounded by an accretion disk. Viscosity within the disk heats it, resulting in thermal 
emission. X-rays are either created in the disk, or in a hot corona above the disk. UV and 
X-ray photons radiate outwards isotropically from the central source. An obscuring torus 
or disk in the plane of the accretion disk absorbs some of the ionizing radiation, resulting 
in an anisotropic radiation field. 

Broad emission lines arise in the region near the SMBH, within the obscuring torus. 
Gas clouds orbiting in the deep potential well are moving quickly, and the lines are doppler 
broadened to ~ 1000 km s _1 . This region is called the "Broad Line Region" (BLR). 
Observers looking through the opening in the torus can see the broad line emission from 
the gas. However, the BLR is obscured for viewers along the plane of the torus. In some 
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Figure 2.2: Schematic geometric model for the unification of radio loud AGN, from 
Urrv Padovanil (jl995l ). The SMBH and accretion disk are at the center. The Broad Line 



Region (BLR) is inside the obscuring torus. The Narrow Line Region (NLR) is exterior 
to the torus and hence visible at all orientations. The radio jets are also shown, forming 
in the vicinity of the SMBH and propagating outwards. Figure Copyright Astronomical 
Society of the Pacific, reprinted by permission of the authors. 

cases, broad emission lines are seen only in polarized light. This is interpreted as due to 
scattering of the broad emission lines off the obscuring torus into the line of sight, implying 
these objects are at some intermediate orientation. 

Narrow emission lines originate from gas clouds in the "Narrow Line Region" (NLR). 
This region of lower density gas is typically farther from the SMBH and so moves more 
slowly, resulting in narrower emission lines. The NLR exists primarily outside the obscur- 
ing torus, and so is not strongly orientation dependent. Forbidden radiative transitions 
occur here. Despite the low values of the Einstein A u \ coefficients, radiative de-excitations 
dominate over collisional de-excitations due to the low density. 

The specific nature of the obscuring geometry is disputed. Initial suggestions were of 
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a atomic or molecular torus (as illustrated in Figure 12. 2p . However, a uniformly dense 
obscuring medium was unable to reproduce observations. Modifications to this mode l 
resulted in a "clumpy" torus formed from discrete clouds (e.g. 



Nenkova et al 



2002 



2008) 



In modeling the observations of Cygnus A (discussed in Chapter [3]) , a model for a clumpy 
torus was tested. An alternate pr oposed obsc uring geometry arises from a wind of material 



blown off the accretion disk (e.g. 



Elvis 



2000) 



The radio structure can also be influenced by the orientation. The jets, moving at rela- 
tivists speeds, are doppler boosted, causing them to appear brighter or fainter depending 
on their relative motion. Evidence for this is seen in VLBI studies of radio cores. When 
both jets are detected, one jet (the "jet") is observed to be brighter than the other (the 
"counterjet"). With the assumption that the jet and counterjet are intrinsically the same 
luminosity luminosity, the inclination of the AGN can be estimated. 

The extended radio emission results from the interaction of the jet with the diffuse 
medium outside the galaxy. The lifetime of relati vistic electrons in the ho tspots is shorter 
than the light travel time to the radio core (e.g. 



Hargrave Sz Rvle 



19741 ). The electrons 



might be re-accelerated in the hotspo ts by 1st order Fermi acceleration before back- flowing 



to create the diffuse lobes seen (e.g. 



Scheuer 



19741 ). Due to the isotropic nature of these 



phenomenon, these features are not likely to be relativistically beamed. This can be seen 
in sources where the hotspot on the side corresponding to the counterjet is brighter than 
the jet side hotspot. 

In this unification scheme, narrow line radio galaxies (NLRGs) are sources viewed 
through the plane of the torus. Radio- loud quasars are those sources viewed nearly down 
the jet axis (ie - through the "hole" of the torus). Broad line radio galaxies (BLRGs) are 
those viewed at intermediate angles, where the broad line region is visible, but the jet is 
not yet significantly doppler boosted. 



14 



2.2 Cygnus A 

Cygnus A is a rare nearby, powerful radio galaxy. At z = 0.05qj it is the most luminous 
radio galaxy out to z ~ 1. Using a cosmol ogy with Hg = 73 km s _1 Mpc -1 , this puts 



20071 ) . The host galaxy is quite 



Cygnus A at a distance of 233.4 ± 16.3 Mpc (jSpergel et al. 
bright at My = —22.6. It is believed to contain a hidden Type I quasar which is heavily 
obscured by dust. Inferences from X-ray observations (fur ther discussed below) suggest 



the AGN's bolometric luminosity is 0.5 — 2 x 10 46 erg s 1 ( Tadhunter et al 



20031 ). 



2.2.1 Observational Characteristics 

Due to the significant radio emission, Cygnus A is classified as a "radio galaxy" and its 
intrinsic radio power is comparable to radio-loud AGN visible in the distant Universe 
(z ~ 1, Dl ~ 6.6 Gpc). As a result, it provides an excellent opportunity to study the 
activity of a powerful AGN and its influence on the host galaxy. Its proximity permits the 
detailed study of the AGN on spatial scales unachievable for its more distant counterparts. 

Cygnus A has been well studied across the electromagnetic spectrum (figure [273]) . This 
section summarizes the observational results, focusing on the nuclear regions, where AGN 
and star formation activity dominate the observed radiation. For a more comprehensive 
rev iew of Cygnus A, which also describes the host galaxy and the extended radio source 



sec 



Carilli fc Barthell fll996|). 
Combining data from the various observations and analyses discussed below, the spec- 
tral energy distribution (SED) of the nucleus of Cygnus A can be constructed. Figure 
shows the flux density at the variety of frequencies for which data are available. 



3 z is the redshift of the source, defined as: 1 + z = A is the observed wavelength of an observed 
(emission) line and Ao is the rest wavelength. 
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Figure 2.3: Broadband spectral energy distribution for Cygnus A. The data were taken 
from a variety of sources in the literature. See text for references. 



Radio Properties 



Cygnus A has be en observed with a variety o f radio observatories at 



frequ encies (e.g. 



Wright &: Birkinshaw 



1984; 



Alexander et al 



many resolutions anc 



1984; 



Meisenheimer et al. 



200 ll ). The extended radio source is a prototypical, edge-brightened FR II (Fig 12.11 right) 
with both the jet and counter-jet visible. It displays obvious hotspots and radio lobes 
similar to those observed in many other triple (core + 2 lobes) radio sources. 

Core radio fluxe s have been obtained from the literature from a few GHz through 



the mm (figure [2 



1984 



Eales et al 



Carilli fc Barthel 



1989; 



Salter et al 



1996; 



Alexander et al 



1984 



Wright Birkinshaw 



19891 ) . At VLA resolutions the unresolved radio core 
has a flat spectrum with spectral index a=0.18 {f u oc f~ a ). However, in high resolution 
VLBI observations, the core spectrum is inverted between 5 and 43 GHz, possibly due 



either to synchrotron self-absorption (SSA) or free-free absorption from a molecular torus 
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(jCarilli et al.lll994l ). The VLA "core" fluxes likely include emission from both the VLBI 
core and the inner jet. 

Both the jet and counter-jet are detected in VLBI observations, with the counter-jet 
being fainter. Assuming the jet and counter-jet are intrinsically the same luminosity (e.g. 



Blandford fc Rees 



19741 ). this indicates the radio source is oriented at some intermediate 



angle to the line of sight. Application of the Symmetric Twin Relativisti c Jet Model to 



these data constrains the orientation of the radio jets to 50° < i < 85° (jSorathia et al. 



19961 ). Following the unification scheme for radio- loud AGN, this suggests that we are 
looking in the plane of the obscuring torus. This is consistent with a geometry in which 
the central SMBH and BLR is hidden from view. 



VLBI observations detected broad HI absorption against the nucleus (jConway &: Blanco 
19951 ). Corresponding to a column density of N# = 2.54 ± 0.44 x 10 19 T sp i n cm -2 (where 



T S pi n is the spin temperature of hydrogen), the absorption has been attributed to the cir- 
cumnuclear torus. The torus is thought to either be atomic (with T~ 8000 K) or molecular 
(T~ 1000 K). 



Conway fc Blancol ()1995l ). the 



Although molecular absorption was not detected by 
observations do not rule out a molecular torus. Radiative excitation of molecules by the 
AGN can result in the absence of molecular absorption lines in lower energy transitions. 



Millimeter / Sub-millimeter 

Millimeter and sub- millimeter observations at two epochs using the S ubmillimetre Commo n 



Robsonetal 



4l998j). 



User Bolometer Array (SCUBA) at the JCMT were presented by 
They find a powerlaw core (likely due to synchrotron radiation, a = 0.6) with a steeper 
powerlaw than observed at lower frequencies. Their data suggests a break or cutoff in the 
synchrotron spectrum at the highest frequencies. The cause may either be depletion of 
energetic electrons, or continuous injection of a powerlaw distribution of electrons, both 
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of which would cause the spectrum to steepen. 

An upper limit was established at 350 fxm, consistent with thermal emission from dust. 
However, the 450 fim flux indicates there is no significant amount of cold dust present. 



Based on these fluxes, 



Robson et al. 



(1998) conclude the emission is dominated by warm 



dust with 37 K < T < 80 K. 



Optical Properties 



Firs t identified by 



cal (Matthews et al 



Baade &: Minkowskil (|1954l ) , the host galaxy of Cygnus A is a cD ellipti 



19641 ) which is also the brightest member of a cluster (see figure 



E3E 



for an optical image of Cygnus A). Ground based observations show the optical nucleus 
has a double structure which is interpreted to be a single nucleus with an obscuring dust 



lane ( Qsterbrock 



1983; 



Tadhunter et al 



19941 ). 



Figure 2.4: Digital Sky Survey (DSS) red image of Cygnus A (center of the frame). 



Tadhunter et al 



(|l994T ) obtained high resolution optical spectra of the nuclear region 
of Cygnus A. Their data were consistent with previous results noting very strong emission 
lines but a relatively weak continuum. The spectra indicate AGN photoionization is the 
dominant ionization source, although they conclude shocks may also contribute. However, 



4 The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Gov- 
ernment grant NAG W-2166. The images of these surveys are based on photographic data obtained using 
the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were 
processed into the present compressed digital form with the permission of these institutions. 
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the coexistence of high and low ionization lines indicates some modifications to AGN pho- 
toionization models are necessary. The ionization field is anisotropic (conical morphology) 
and aligned with the radio axis, consistent with unified models for radio- loud AGN (see 
£ j2.1.ip . Also, in order to explain the observed nitrogen ratios, abundances must be much 
higher than solar values. 

Cygnus A i s thought to contain an ob scured quasar similar to those seen at high 



redshift. While 



Jackson &; Tadhunterl (|1993l ) were unab 



light, it was detected in later observations by 



Ogle et al. 



e to d etect broad Ha in polarized 



(|l997l ). The polarized broad lines 



were observed within the ionization cone noted above, with a possible scattering origin 
in the NLR. The detection of broad Ha emission lends support for the existence of an 
obscured BLR. 



Tadhunter et al 



((2003) obtained Hubble Space Telescope (HST) Space Telescope Imag- 
ing Spectrograph (STIS) spectra of the nuclear region. Stepping the slit across the nucleus, 
a velocity gradient indicative of rotation around the radio axis was observed. Modeling 
the velocity as due to the potential of a SMBH and stellar mass distribution (measured 
from a 1.6/iin NICMOS image) gives a SMBH mass of 2.5 ± 0.7 x 10 9 M^. This measure- 



ment is consistent with b 



Ferrarese fe Merritt 



200ol : 



ack hole mass - host galaxy relations (jMagorrian et al 



Gebhardt et al 



2000) 



1998 



HST imaging has also revealed star fo rmation in the central region of Cygnus A, which 



began < 1 Gyr ago (jJackson et al 



1998). Lying in a 4kpc ring around the nucleus, that 
is oriented orthogonal to the radio axis, emission from the starburst is detected in the 
Spitzer observations described below. 



Infrared Properties 



Keck NIR imaging with adaptive optics was presented by 



Canalizo et al 



(120031 ). They 



find a secondary point source 0.4" from the nucleus (K~ 19). The IR-Optical SED of 
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this off-nucleus point source is consistent with a smaller galaxy which has been tidally 
stripped of much of its original mass. It is likely a dense core of older stars, which is 
all that remains from the merging galaxy. The merging scenario can explain the recent 
(~ 10 7 yr) triggering of AGN activity, the disturbed kinematics seen in Cygnus A, and 
the dust lane across the nucleus. 



Tadhunter et al 



(11999 ) obtained NICMOS images of Cygnus A from HST. A point 



source in their image 



Tadhunter et al 



(12003). 



Keck K-band spectra of the nuclear region were obtained by 
While they did not provide a full analysis, spatially extended Paa and H2 emission was 
detected, and fits to these lines were used to constrain the SMBH mass (Af, = (2.5 ± 
0.7 ) x 10 9 Mg). 



Bellamy fc Tadhunter! (|2004l ) presented a more detailed analysis of the K-band spectra, 
finding complicated emission line properties suggesting an infalling molecular cloud. The 
location and velocity of the molecular (#2) emission is inconsistent with interpretation as 
an outflowing cloud. The H2 line ratios are consistent with X-ray heating, indicating the 
cloud is likely excited directly by the X-rays, via hard X-rays which can pe netrate the ob- 



scurin g: torus. The interpretation as an infalling cloud is consistent with the 



Canalizo et al. 



(2003) picture where Cygnus A is currently undergoing a minor merger. 

Additionally, the H2 line was seen in several components, both red- and blue-shifted 
relative to the systemic velocity. This was interpreted as emission from a rotating torus. 
The observed line ratios are consistent with excitation by X-rays (from the AGN), while 
likely ruling out shocks as a possible excitation method. 

The wide range of flux values observed in the near- IR (figure E3 ~ 10 14 Hz) illustrates 
the importance of aperture matching when constructing such a broadband SED. The 
2MASS aperture encompasses much of the host galaxy, including much stellar emission. 
The Keck K-band aperture is smaller, while the NICMOS observation is the point source 
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continuum flux associated with the nucleus. 

Mid-infrared (11.7 fim) imaging of Cygnus A were obtained by 



(200 



Whysong &; Antonucci 



Based on the observed nuclear flux, and an assumed spectral shape similar to 



those for radio-loud Palomar-Green quasars, they estimate the bolometric luminosity of 
the AGN in Cygnus A to be ~ 1.5 x 10 45 erg s _1 . 

UV Observations 



Zirbel fc Bauml p998 ). Undetected 



HST UV observations of Cygnus A were presented by 

at the shortest wavelengths observed (F130IvJf , F152M, F170M, & F220M), the F320M, 
F342W, and F372M data show no nuclear component in modeling, consistent with the 
obscuration seen in the optical. Any UV flux from the AGN would be absorbed by the 
optically thick torus and re-radiated in the infrared. 



Using the Faint Object Camera on HST, 



Antonucci et al 



fll99J ) found broad (FWHM 



1 o 

~ 7500 km s ) [Mg II] 2799A emission, supportive of the notion that Cygnus A contains 
a hidden quasar. 



X-ray Properties 



Cygn us A has b een observed with a variet y of X-ray i nstruments, includ i ng Uhuru 

2s, 



1972), ROSAT (Reynolds & Fabian 



and INTEGRAL ( 



Beckmann et al 



1996), Chandra (lYoung et al 



200^ 



Gur sky et al 



Evans et al 



2006) 



2006). These observations support the existence of the 



hidden quasar in Cygnus A based on the appearance of an unresolved hard X-ray spectrum 
coincident with the radio core. 



While of relatively poor angular resolution, ROSAT observations by ([Reynolds Fabian 



19961 ) were able to resolve the hot ICM around Cygnus A. Spectral fitting gives tempera- 



5 Hubble Space Telescope filters are designated by F, the central wavelength (in nm), and 'N' for a 
narrow filter, 'M' for a medium width filter, and 'W for a wide filter. 
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tures of ~ 3 — 4 keV over the cluster. Also, evidence for a cooling flow of ~ 250 M Q yr 1 
was seen. 



Young et al 



(2002) conducted Chandra ACIS observations of 32ks and lOks. The 
X-ray flux above 2keV was peaked at the nucleus while lower energy flux was more dis- 
tributed. The hard X-ray nucleus is coincident with the radio nucleus and is less than 0.4" 
in size. Modeling the emission gives a neutral column density Njj = 2]to' 2 x 10 23 cm~ 2 
towards the nucleus. 

I NTEGRAL observatio ns from 20-f00 keV give a luminosity of log L 2 



s 1 (Be ckmann et al 



44.71 erg 



2006) integrated over the galaxy. The resolution of the imager is 
~ 12 arcmin, and so this value may be subject to contamination from the hot X-ray gas 
associated with the cluster. 



2.2.2 Summary: Current Knowledge 

Based on the abundance of multi-wavelength data, a general picture can be constructed 
for Cygnus A. It is a large elliptical galaxy in a poor cluster. The galaxy hosts a SMBH 
with mass typical of SMBHs in other galaxies with similar large-scale properties. It is 
currently undergoing a minor merger with a smaller galaxy which it has tidally stripped. 
The existence of a dust lane and infalling gas can be connected to this merger event 
which is likely linked to the nuclear activity (AGN and starburst). Hard X-ray emission 
coincident with the nucleus is likely due to emission from the accretion disk surrounding 
the SMBH. 

While the specifics of physical components derived from modeling will be discussed in 
Chapter O it is useful to briefly mention the major features of the SEP 12.31 At the lower 
frequencies (radio through sub- mm, 10 9 — 10 12 Hz), the SED clearly follows a powerlaw 
- indicative of non-thermal synchrotron emission from the radio core. At sufficiently low 
frequencies, the radiating plasma is expected to become optically thick and synchrotron 
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self-absorption will occur. 

The powerlaw is interrupted in the FIR and thermal emission is seen. This ther- 



ma 



bump, which extend s almost to the NIR, has an effective temperature of ~ 75 



K (|Diorgovski et al 



199ll ). Emission from the circumnuclear starburst and/or torus- 



reprocessed UV and X-ray radiation from the AGN contributes to this bump. At ~ 10 /im, 
the extrapolated powerlaw would exceed the emission seen, so the synchrotron spectrum 
should break somewhere in the IR, likely be due to a depletion of the relativistic electron 
population at higher energies. 
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Chapter 3 



Observations 



Infrared observations of AGN are critical in understanding the nature of the emission. 
In the case of Cygnus A, observations in the mid and far infrared can shed light on the 
relative importance of the AGN activity and star formation. These phenomena have 
been studied with previous space-based infrared missions including the Infrared Space 
Observatory (ISO) and the Infrared Astronomy Satellite (IRAS). 

Technological improvements since these missions enable a more detailed study of the 
nuclear regions of Cygnus A. The following pages discuss the Spitzer Space Telescope and 
the spectroscopic observations utilized in constructing and modeling the spectral energy 
distribution of Cygnus A. 



3.1 Spitzer Space Telescope 



Launched in 2003, the Spitzer Space telescope i mproves on previous 



space-based infrared 



2004). While its primary 



observatories in both resolution and sensitivity (j Werner et al. 
mirror is only slightly larger than previous space-borne infrared missions (at 85cm), Spitzer 
features greatly improved detectors. Additionally, the combination of passive and active 
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cooling should drastically improve the usable lifetime of the telescope. 
There are three instruments on boar d: the Infrared Arr ay Camera (IRAC. 



Fazio et al. 



20041 ). the Infrared 



Photometer (MIPS. 



Spectrograph (IR S, iHouck et al 



Riekeetal 



2004 ) . and the Multiband Imaging 



20041 ). These instruments collectively cover a large portion 



of wavelengths between 3.6 and 160 fj,m in both imaging and spectroscopy modes. 

The optical design is a cassegrain telescope, utilizing beryllium optics. The mirror is 
cooled to ~ 5.5K. This cooling is achieved through the separation of the satellite into a 
"warm assembly" and a "cold assembly". The warm assembly contains the communica- 
tions equipment and warm electronics for the instruments. The sun shield is also attached 
here. 

The cold assembly is attached to the rest of the satellite through low thermal conduc- 
tivity struts and insulated with Mylar thermal blankets and a radiation shield. To conserve 
cyrogens a "warm launch" was utilized in which most of the satellite was launched at am- 
bient temperature. Once in the Earth-trailing orbit, the telescope was allowed to passively 
cool, after which the cyrogens were used to further cool the instruments to their operating 
temperature. A satellite orbiting Earth is exposed to both reflected sunlight as well as 
thermal radiation from the warm Earth. The Earth-trailing orbit allows the satellite to 
avoid the thermal consequences of close proximity to the Earth, conserving cryogens and 
extending the mission. 

The focal plane of Spitzer is shown in Figure I3TT1 (not to scale). The slit orientations 
of both the SL and LL modules are shown, as well as the locations of the various imaging 
arrays. 
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Figur e 3.1: A d iagram (not to scale) of the focal plane of the Spitzer Space Telescope. 
FromlSOMj (|2006l ). 



3.2 The Infrared Spectrograph 



The Infrared Spectrograph (IRS) on Spitzer is primarily designed to obtain spectra from 
~5-38^m. It is capable of medium- (R ~ 600) and low-resolution (R ~ 60-120) spectro- 



scopic observations as well as "peak-up" imaging (jHouck et al 



200J). Here I will only 
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be concerned with the properties of the low-resolution observing mode. Details of the 
medium-resolution and peak-up imaging modes can be found in the references. 

IRS has two modules covering the low-resolution mode. The Short-Low (SL) module 
observes the short-wavelength (~5-14 /im) region while the Long-Low (LL) module ob- 
serves the long wavelength region (~14-35 /im). There is overlap between the SL and LL 
modules for consistent calibration between the two spectral regions. 

The SL slit is 3.7" wide with a corresponding pixel size of 1.8 " pixel -1 . The LL slit is 
oriented at an angle of 96.5° to the SL slit. It is 10.7 " across with a pixel scale of 5.1 " 
pixel -1 . 

Each module covers three spectra orders. Table 13.11 lists the wavelength coverage of 
each order. In both the SL and LL modules, the third order is a "bonus" order, included 
to provide additional data at the overlap between the two primary orders. 

Table 3.1: Wavelength Coverage of the IRS Orders 



Order A (/im) 



SL2 


5.2 


- 7.7 


SL1 


7.4 - 


- 14.5 


SL3 ("bonus") 


7.3 


- 8.7 


LL2 


14.0 


- 21.3 


LL1 


19.5 


- 38.0 


LL3 ("bonus") 


19.4 


- 21.7 



3.2.1 Detectors 

The IRS detector is silicon, doped with arsenide (As) for short wavelengths, and antimony 
(Sb) for the long wavelengths. The Si:As chips are sensitive from 5-26 fim and the Si:Sb 
chips from 14-40 /im. The detective quantum efficiency is a function of both wavelength 
and the bias voltage. For a V^ as =1.5V the DQE for the Si:Sb chips is ~15% (ranging 
from 5-25% over 20-36/im). 

In order to avoid saturation of the detector and to lessen the impact of cosmic rays 
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Figure 3.2: Layout of the detector for the low resolution IRS modules. LEFT: SL module. 
The strip to the left is the SL1 order. The middle section includes the SL2 and SL3 orders. 
To the right are the two peak-up imaging arrays. Right: LL module. 

(CRs) the detectors are read out in a ramp fashion. The total integration is split up into 
discrete sample times defined by the user. The pixel values are read out at the end of each 
sample time and recorded (Figure I3.3[) . Non-destructive readout of the detector enables 
the voltages to be recorded multiple times. 

If a pixel becomes saturated during an observation (either from a source, or a CR 
hit), the voltage will remain at a (roughly) constant level for the rest of the observation. 
Measurements of the ramp slope until saturation enable a determination of the flux. 

Cosmic rays hits that don't saturate the pixel will result in significant drops in voltage, 
and the slopes on either side of the impact can be used to determine the actual incident 
flux. 
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Figure 3.3: Figure demonstrating the data collection method of IRS. See text for details 
(jSOMll2006h . 



3.2.2 IRS Observations 

IRS provides two spectroscopic observing modes: "staring mode" and "mapping mode". 
In staring mode, the slit can be positioned directly on the source for a single observation. 
In "mapping mode", the slit is stepped across the source in increments equal to half the 
slit width. This Nyquist sampling enables reconstruction of the spatial distribution of 
the emission. The Spitzer Infrared Nearby Galaxy Survey (SINGS) Legacy program has 
written and provided software to take IRS spectral mapping data and convert into a data 
cube with both spatial dimensions on the sky (RA & Dec) as well as a wavelength axis. 
The "Cube Builder for IRS Spectra Maps" (CUBISM) software can use IRS mapping 
observations to rec onstruct the brightness distribution of sources across a ~30 /um range 



( Smith et al. 



20071 ). In addition, images can be constructed in specific emission lines, 
spatially localizing the emitting region. 

In order to spatially map the MIR emission from Cygnus A, observations were per- 



formed in mapping mode. 13 observations were made with the SL module, stepping across 
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the source in 1.8" increments (50% of the slit width). Similarly, five observations were 
made with the LL module, again stepping across the source in half-slit width increments 
(4.85"). In both cases, off source data for background subtraction were obtained roughly 
100" perpendicular to the slit stepping direction. These off-source data were obtained 
automatically as part of the observing program. During an observation utilizing the LL 
slit on source, off-source observations are simultaneously made with the SL slit, and vice 
versa. Figure 13.41 shows the slit locations across the source. 




Figure 3.4: 2MASS K-band image of Cygnus A with IRS slit positions overlaid. The 
narrow boxes represent the SL slit, and the wide boxes show the LL slit positions. Note 
that the slits were stepped across the source in half-slit-width increments. The off-source 
LL slit positions are not shown, but are to the south-east. 
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3.3 Data Reduction 



Data from the Spitzer Space Telescope is processed in two stages. First, a software pipeline 
at the SSC performs standard image processing. These Basic Calibrated Data (BCD) 
products are then made available to observers for further reduction and analysis through 
the Spitzer Science Center (SSC). The IRS instrument team and the SSC have provided 
software tools for final end-user processing. 

3.3.1 Pipeline Reduction 

After being downlinked from the satellite individual observations are subject pipeline 



reduction to remove most instrumental effects and CR hits from the data. 



will p rovide 



2006) 



a short overview, however for more detail, see the IRS Data Handbook (IRS ! 

The first step in the pipeline flags pixels affected by CRs. This is done by locating 
jumps in ramp readouts. Due to the non-destructive readout, the ramp slope can be 
measured on either side of the jump, and the correct incident flux can be recovered (see 
Figure 13.31 and text above for more details) . 

On these detectors the charge on a pixel is affected by the overall flux on the array 
(called "droop"). The next step in the pipeline factors in the total flux incident on the 
array and corrects individual pixel values. Next, dark frames are subtracted from the 
data, to remove charge due to thermal excitation. As the pixel response to incident flux is 
generally non-linear, the pixel responses are linearized using a model for each individual 
pixel. 

With cosmic ray hits removed and the pixel response linearized, the incident flux can 
be determined by measuring the slope of the ramp. 

The IRS chips feature 4 readout channels and slightly varying gains in each channel 
can cause a "jailbar" effect in the data. This is corrected next in the pipeline. The final 
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corrections involve removing stray light in the SL regime and applying a flat field to correct 
for overall pixel sensitivity variations. 

The output of the pipeline is a set of files comprising the basic calibrated data (BCD). 
They are downloaded by the end user and are the starting point for further calibration 
and analysis. The primary BCD files are: 

1. OBSJDENTIFIER±> cd . fit s 

2. OBSJDENTIFIERJuncMs 

3. OBSJDENTIFIER_bmaskMs 

The OBS_IDENTIFIER.bcd.tits file contains the pipeline processed observations. It 
is laid out as shown in figure [3T2l OBS-IDENTIFIER _func.fits contains estimates of the 
uncertainty in the pixel values. Finally, OBS-IDENTIFIERJomask.tits is a mask file of 
pixels (flagged both by the pipeline and the end user). 

The pipeline also provides post-BCD data products in which spectra have already been 
pipeline extracted and flux calibrated. However, the post-BCD products have limited 
flagging of bad pixels and are generally only suitable for a first look at the observations. 
They can also be used to examine the data for potential issues and plan processing of the 
BCDs. 

In addition the above processed products, the SSC provides the raw data files (pre- 
pipeline reduction) to enable end-users to process the data in a non-standard fashion. 

3.3.2 BCD Reduction 

Prior to extracting and calibrating the spectra, additional instrumental effects not ad- 
dressed by the pipeline processing must be removed. The detector array contains pixels 
which are defective or have non- linear responses (even after pipeline calibrations). An IDL 
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script (IRSCLEAN) is available to assist in eliminating bad pixels from the BCD files. It 
utilizes both bad pixel masks from the SSC and user flagged pixels to build a mask file 
used by extraction routines to ignore bad pixels. 

Two software products are available for calibration and extract ion of IRS spectra. The 



SSC provides the "Spitzer IRS Custom Extraction" (SPICE) tool (ISPI 



2006h . and the IRS 



instrument team at Cornell has developed a set of tools called the "Spectroscopy Modeling 



2004). Both tools can be used to 



Analysis and Reduction Tool" (SMART) (jHigdon et al 
extract and flux calibrate IRS spectra. 

Once the BCD files have been "cleaned" (using IRSCLEAN), the background must 
be subtracted from the images. For the low-resolution modules, this can be accomplished 
using off-source data. When a source with SL1, the SL2 slit is off source by ~ 100". 
By subtracting the SL2 off source FITS file from the SL1 on-source FITS file, the local 
background can be subtracted, leaving just the spectra. 

After the background has been subtracted, the data was loaded into SPICE for ex- 
traction and calibration. Three files are required for the extraction: _bcd.fits, _func.fits, 
and _bmask.fits. The first file contains the observation, the second has uncertainty infor- 
mation, and the final file contains information about bad and "rogue" pixels (output from 
IRSCLEAN). "Rogue" pixels are those pixels which are not consistently bad, but have 
anomalous values in a particular exposure. The first step in SPICE is to determine the 
location of the spectra. In the case of a point source, a profile in the position direction 
will show a peak at the location of a source. The extraction aperture is centered on this 
peak. The aperture scales with wavelength to account for the changing size of the point 
spread function (PSF). 

For a point source, the aperture is wavelength-dependent, due to the wavelength de- 
pendent size of the PSF. The pixels are summed across this aperture to give counts vs 
wavelength. An additional wavelength dependent correction is applied to account for light 
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that is scattered out of the slit (slit loss correction). 

This slit loss correction is relatively straight-forward for a point source, where the 
correction for light not incident on the slit can easily be made. However, in the case of 
an extended object, flux can also be scattered into the slit. Any correction must include 
both light scattered out of the slit due to the PSF, as well as light scattered into the slit 
from the PSF of neighboring source regions. This can only be properly done if the spatial 
distribution of the source is known. SPICE includes an extended source extraction which 
assumes an extended source of uniform surface brightness. In this case, as much light is 
scattered into the slit as out, and a "slit loss correction function" (SLCF) uncorrects the 
assumed amount of light scattered into the slit in the case of a point source. 

Finally, the BCD values are flux calibrated in units of Janskys (1 Jy = 10 -23 erg 
s _1 cm~ 2 Hz -1 ). The flux calibration is obtained through observations of standard stars. 
Stellar atmosphere s have been modeled for these stars and facilitate the detailed calibration 



(Decinetal 



2004]). 

This process much be repeated for both orders of each module, resulting in multiple 
output files which must be stitched together for a final spectrum. The modules overlap 
in wavelength, so they can be matched in flux (for a point source). The "bonus" order 
in both modules is extracted with the I s4 order and can be used to provide an additional 
constraint. 

As noted earlier, the SL and LL modules have markedly different slit sizes. For point 
sources, the amount of light scattered out of both slits can be modeled and calibrated. 
In the case of an extended source, the wider LL slit samples a much larger area on the 
source. This difference in spatial sampling is evident in a significant jump in flux seen at 
the intersection of the SL and LL modules (see Fig [33]). 

Spectral mapping observations offer a solution to this problem. By stepping the slit 
across the source (at Nyquist limit), a data cube can be constructed (two spatial axes, 
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Figure 3.5: Demonstation of the flux discontinuity near 14/xm in SL and LL observations of 
an extended source (Cygnus A). Note that this jump is present regardless of the extraction 
method used (point source or extended source). 

and one wavelength axis). A spectrum can then be extracted from the cube of any arbi- 
trary aperture at any location. This facilitates the extraction of spectra across the entire 
wavelength range using the same size aperture. 

3.4 Cygnus A IRS Data Reduction 

Cygnus A was observed using both the SL and LL modules in mapping mode. The slit was 
stepped across the source in half-slit width increments to facilitate the construction of a 
data cube. SL observations were taken at 13 slit locations (1.8" shifts) and LL observations 
were taken at 5 slit locations (4.85" shifts). For all four modules, the effective integration 
time was 14s at each slit location. 

After pipeline processing at the SSC, the individual BCDs were cleaned of bad and 
warm pixels using IRSCLEAN. The cleaned BCDs were then loaded into CUBISM where 
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Table 3.2: Additional Infrared Flux Densities from Spitzer 
A ({j,m) F u (Jy) 



4.5 0.010 ±0.003 

8 0.054 ±0.013 
70 2.20 ±0.11 

160 0.668 ± 0.033 



they were background subtracted and assembled into cubes for the SL and LL modules. 
A 20" circular aperture was used to extract the IRS spectrum from each cube (shown in 
figure [3T6]) . A 20" aperture was chosen to cover a well sampled region of sky. The LL 
slit has a 10.7" width, so the aperture was chosen to be as small as possible while still 
properly sampling the reconstrctued data cube. 

3.5 Additional Spitzer Data 

Imaging was also done using Spitzer's Infrared Array Camera (IRAC) at 4.5 and 8.0/mi (PI: 
Harris). The unpublished data were downloaded from the archive, reduced and calibrated 
according to the instrument manual and included on SED and in the modeling. The flux 
de nsities are given i n table 13.21 



Shi et al 



(2005) observed Cygnus A using the Multiband Imaging Photometer for 
SIRTF (MIPS) instrument on Spitzer at 24, 70, and 160/im. Their 24/im flux is consistent 
with our IRS observations. However, their measurements between 70 and 160^mi give a 
decrease in flux steeper than the Rayleigh-Jeans tail. As the flux in the MIR is probably 
due to thermal dust emission, such a steep falloff is likely unphysical. The BCDs were 
downloaded from the SSC and re-reduced. The new reduction resulted in a higher flux 
density measurement at 160 /um, more consistent with thermal emission. The new flux 
densities obtained from the archival MIPS data are given in table I3T21 
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Figure 3.6: Top: IRS spectrum for Cygnus A, as extracted from a 20" aperture centered 
on the nucleus. Bottom: Same as top, but on a log-linear scale. 
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Chapter 4 

Modeling Approach 



As is evident from §2.2.2l the spectral energy distribution of Cygnus A is quite complicated, 
featuring both thermal and non-thermal emission from a variety of physical regions on 
drastically different size scales (ranging from sub-pc to kpc scales). The combination 
of the various physical processes at work (AGN activity, starburst) creates a possibly 
complicated infrared spectrum. Here only the continuum emission will be modeled. 

Strong emission lines are clearly evident in the IRS spectrum. With the exception of the 
polycyclic aromatic hydrocarbon (PAH) features, these emission lines generally result from 
either photoionization from the AGN or shock-ionization resulting from inte raction with 



the r adio jet. The emission lines and PAH features were fit using PAHFIT (jSmith et al 



20071 ). Using the fits, the emission lines were removed from the spectrum, leaving only the 
underlying continuum and PAH features. 

The resulting spectrum of continuum + PAH features was fit using a combination of 
AGN and Starburst processes as described in ^4.11 This choice of models was motivated 
by the observed properties previously described (Chapter [2]). 

The cleaned IRS spectrum was combined with a variety of data from the literature 
to produce a spectral energy distribution from the radio through the near IR. All fluxes 
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Figure 4.1: Cygnus A IRS spectum with the emission lines removed. 



were measured at the nucleus. Where possible, fluxes were extracted from a 20" aper- 



ture 


to 


match 


that 


(1989 


), 


Salter et al. 


Robson et al. 


( 


1998) 



(1989) 



RS observations. Data were obtained from: 



Alexander et al. 



(|1984l ) 



Wright Sz Birkinshawl ( 



Eales et al. 



1984 ). and 



19981 ). Table 14.11 lists the values, aperture sizes, and references for the 



data points. 



Table 4.1: Cygnus A fluxes adopted from the literature. 

v (GHz) /„ (Jy) Ap (") Ref. 

2.7 1.5 ±0.4 3.7x5.8 Alexander et al. (1984) 

5 0.95 ±0.2 2.0x3.1 Alexander et al. (1984) 

15 1.22 ±0.2 0.7x1.0 Alexander et al. (1984) 

89 0.7 ±0.07 1.5x2.8 Wright fc Birkinshaw (1984) 

150 1.11 ±0.1 33.5 Robson et al. (1998) 

220 0.65 ±0.1 22.9 Eobson et al. (1998) 

230 0.59 ±0.07 11 Salter et al. (1989) 

272 0.58 ±0.06 19 Eales et al. (1989 ) 

350 0.46 ±0.06 14.4 Eobson et al. (1998) 
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4.1 Model Choices 



Initially, a simplistic model consisting of blackbody emission with an intervening screen 
of dust was fit to the data, with very poor results. A two temperature component model 
was also fit, with similarly poor results. Based on this, I formulated more detailed and 
specific models. 

The choices of models have been dictated by evidence from other observations of the 
nuclear activity in Cygnus A. The modeling is aimed at determining the relative contri- 
butions from each component in the nucleus. Additionally, an estimate of the bolometric 
luminosity of the AGN can be obtained, testing the hypothesis that Cygnus A harbors an 
obscured quasar. 

Thus for the purposes of modeling the SED of Cygnus A, three components have been 
included: 

1. Clumpy, dusty torus 

2. Starburst 

3. Synchrotron emission from a jet 

Each component will be discussed in turn, describing the motivation for its selection in 
the context of Cygnus A and the parameters available to match observations. 

4.1.1 Torus Model: CLUMPY 

Cygnus A is a radio-loud active galaxy believed to be powered by a supermassive black hole 
at its center (Chapter [2]). Observational evidence indicates the inner region of the AGN 
is obscured, possibly by a geometrical structure similar to that predicted by unification 
models (see § 12.1. ip . Models describing the predicted AGN structure can be compared 
with the observations to test these unification theories. 
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In Cygnus A, the radio jets are oriented roughly in the plane of the sky, so doppler 
(de)boosting effects are roughly equivalent for both the jet and counter-jet. The primary 
test of unification models in Cygnus A is studying the geometric structure obscuring the 
central SMBH and accretion disk. The original concept of an obscuring torus was of a 
thick "donut" surrounding the accretion disk, obscuring the accretion disk from certain 
viewing angles. The stability of such a structure was called into question and simulations 
indicated that a thick molecular torus would collapse into a disk within a few dynamical 
timescales. 

Improvements in the torus models have resulted in a "clumpy" torus, where discrete 
clouds of gas form a torus like object, surrounding t he accretion disk. Ra diation pressure 



ca n support the torus ag ainst gravitational collapse (IPier h Krolik 



Nenkova et al 



1992 ). 



(|2002l ) have constructed a set of clumpy torus models using the DUSTY 
radiative transfer code. Called CLUMPY, their grid of models is available onlinqj. The 
geometry of the model is shown in Figure l-Ol 

to observer 
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Figur e 4.2: Geomet ry of the CLUMPY t orus model as described in lNenkova et al.1 ([2002 
2008). Figure from Nenkova et al. ( 20081 ). used with permission. 



The parameters for the model are given here and explained below: 



: http : //newton.pa.uky . edu/~ clumpy web/ 
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1. N: Average number of clumps along the line of sight 

2. Y = J^: Ratio of outer to inne radius of torus 

3. q: powerlaw index of the radial distribution of clumps 

90— i 

4. a (degrees): characteristic width of the angular distribution of clumps, oc e~ ~ 

5. Ty- Optical depth of a clump at 0.55^/m (all clumps are assumed to have the same 
optical depth) 

The outer radius of the torus R is Y times the inner radius (Rd)- The inner radius is 
determined from a dust sublimation temperature of T = 1500 K (R^ = 0.4*L4j 5 pc, where 
L45 is the bolometric luminosity of the AGN in units of 10 45 erg s _1 ). All models fit used a 
torus with a gaussian angular distribution, with a parameterizing the width of the angular 
distribution from the plane, q is the index of the powerlaw distribution of clumps in the 
radial direction: r~ q . Table B~2l gives the range in values for these parameters during the 
modeling described here. 

In addition to the above parameters the models are provided with and without the 
AGN emission directly visible. There is a finite probability the AGN emission might be 
visible through the clumpy torus. Figure ET31 shows the effect of varying each parameter. 

The AGN spectrum used in the clumpy models is: 



A/a oc < 



A 1 ' 2 


A < X h 


const 


Ah < A < X u 


A -0.5 


A n < A < Xrj 


A- 3 


Xrj < A 



(4.1) 



Where Xrj is the wavelength at which the spectrum becomes Rayleigh-Jeans (1/im 
here). Xh was taken as O.Olfim and X u = 0.1 fim. 
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The inclination (i) determines the degree to which an observer sees the AGN emission 
absorbed by the torus. As % increases, the observer goes from a (mostly) unobscured view 
of the central regions to an edge-on view of the torus. Correspondingly, more optical and 
near-infrared emission is absorbed at higher inclinations. This absorbed radiation is re- 
emitted isotropically, so the far-infrared emission varies little as the inclination is changed. 
At wavelengths shorter than 20 fim , the clouds begin to become optically thick and the 
radiation is no longer isotropic. 

Two files are provided for each combination of parameters: one with only the torus 
emission, and one with the torus emission and powerlaw emission from the accretion disk. 
They are the same at long wavelengths, as the torus dominates. However, in the NIR, the 
accretion disk emission (if directly visible) begins to dominate. 

The characteristic width of the angular distribution (cr) effectively determines the 
"thickness" of the torus in the vertical direction. Increasing this value results in an in- 
creased number of clumps along the line of sight for a given inclination. The additional 
absorption results in additional radiation being absorbed and re-emitted in the far-infrared. 

Varying the powerlaw index of the clump radial distribution (q) effectively changes 
the solid angle subtended by clumps as seen by the AGN. Thus, for a steeply falling 
distribution (large q), the N clumps along a line of sight would tend to be closer to the 
AGN, and thus subtend a larger solid angle, absorbing more radiation. This is manifested 
as a shift in the peak flux towards shorter wavelengths as well as a decrease in flux at 
longer wavelengths with increasing q. Also, silicate absorption around lO^um becomes less 
pronounced. Clouds closer to the AGN are warmer than clouds out, so a more centrally 
peaked distribution results in the overall emission peaking towards shorter wavelengths. 

Varying the size of the torus (Y) has a comparatively smaller effect on the overall 
spectral shape. Small values of Y have slightly higher flux at shorter wavelengths, with 
a corresponding decrease in flux at longer wavelengths, when compared to large values of 
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Y. Large values of Y provide a larger range in individual cloud temperatures, broadening 
the emission. The effect is not very pronounced in figure 14.31 because of the steep radial 
distribution (q = 3) used, resulting in few clouds at large radii. 

Varying the number of clumps along a line of sight (N) changes the overall spectral 
shape. Increasing the number of clumps increases the depth of the silicate absorption as 
well as reprocessing optical and near- infrared radiation into the far- infrared. Increasing 
the optical depth of a clump (-ry) has a similar effect. For larger values of N and ry, the 
silicate absorption becomes saturated and begins to fill back in. 

For the each set of parameter values, the radiation observed at i between 0°-90° is 
computed. 

4.1.2 Starburst Model: Siebenmorgen & Kriigel 



Siebenmorgen Kriigell (|2007l ) have developed a set of starburst models with a variety 
of physical conditions. Assuming spherical symmetry and an ISM with dust properties 
characteristic of the Milky Way, they generated a grid of over 7000 models using radiative 
transfer calculations. The free parameters for the model are: 

1. Starburst radius r 

2. Total Luminosity Lt t 

3. Ratio of bolometric luminosity in O and B stars to total luminosity foB 

4. Visual extinction from the center to the edge of the nucleus A v 

5. Dust density in hotspots around O and B stars n 



Ranges for each parameter are given in table I4~2l 

Emission is broken down into two components: an old stellar population uniformly 
distributed through the volume and hot, luminous O and B stars embedded in dusty hot 
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spots. O and B stars are massive stars (~4-100 M@) with high luminosities (~ 10 2 — 10 6 
Lq) and effective temperatures (T e ff ~2-6x 10 4 K ). As a result they produce a large 
quantity of ionizing photons. 

The density of OB stars is centrally peaked although the model output is the emission 
integrated over the entire starburst. The model grid is available online. 

Figure |4~41 shows the model dependence on the various parameters. The visual extinc- 
tion (Ay) quantifies the radiation absorbed at V-band (550 nm). Increasing Ay results 
in increased absorption at UV to mid-IR wavelengths (A < 70/xm) and increased depth of 
the silicate absorption at 10/im. The energy absorbed at short wavelengths is re-emitted 
in the far-infrared, so increasing Ay results in an increase in far-infrared flux. 

Increasing the dust density in hot spots (n) has a similar effect to increasing Ay, how- 
ever the energy is re-emitted at shorter wavelengths (~ 70//m compared to ~> 100/um). 

The fraction of OB stars has a smaller effect on the spectrum, with an increased foB 
resulting in increased mid-infrared flux and a decrease in optical and far-infrared flux. 

The integrated luminosity of the starburst (Ltot) determines primarily the vertical 
position of each model. There is some effect on the releative strength of the emission 
features as a function of Lf t, with shorter wavelength features becoming stronger at higher 
luminosities. The radius of the starburst (r) has no effect on the spectrum at wavelengths 
shorter than ~ 20/xm. At longer wavelengths, increasing r increases the far-infrared flux. 
This can be attributed to a larger volume of warm dust. 

4.1.3 Synchrotron Emission 

Extrapolating the powerlaw emission in the radio portion of the SED for Cygnus A to the 
mid-IR indicates the powerlaw must break somewhere in mid-IR. The break frequency is 
apparently hidden by the thermal bump, and the contribution of the synchrotron radiation 



^http: //www . eso . org/~rsiebenm/sb_models/ 
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to the sub-mm and IR flux is uncertain. This component is modeled as a self-absorbed 
power-law spectrum which breaks at high frequencies due to a cutoff in the distribution 
of particle energies. At lower frequencies the emitting plasma becomes optically thick 
and results in synchrotron self-absorption (SSA). This turnover at lower frequencies is not 
evident so this portion was not modeled. The accepted functional form is: 

F v oc (l - e-m ai -" 2 ) e -* e -"M (4.2) 

where v Q is the frequency at which the optical depth of the emitting plasma is equal 
to unity, v c is the frequency corresponding to the cutoff in the energy distribution of the 
particles, ol\ is the spectral index in the optically thick regime, 02 is the spectral index 
in the optically thin portion, and r r is the dust screen between the synchrotron emitting 
region and the observer. 

The turnover due to self-absorption of synchrotron radiation occurs at frequencies 
lower than those observed, so «2 and v Q were not free parameters for fitting. In Cygnus A, 
the spectral index a.2 can be measured from core radio fluxes. Additionally, the amplitude 
of the powerlaw is fixed from the same observations. The only unconstrained parameter 
is the break frequency v c which will depend on the contributions from the previous two 
model components to the mid-IR spectrum. 

An alternate model for the synchrotron emission at higher frequencies is a broken 
powerlaw. Aging of the population of relativistic electrons results in a broke n powerlaw 



whos e spectral index steepens at frequencies higher than the break frequencies (jKardashev 



1962j). The functional form adoped is. 



F u oc < 



V 1 V < Vbreak, 

(4.3) 

V~ a2 V > Inbreak 
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In the case of Cygnus A, the spectral index post-break would be 02 = 1.24. 
Fits were run using either the exponential cutoff or the spectral break. Each implies 
fundamentally different populations of relativistic electrons. 

To simplify references to these two fits, they will be referred to as: 

• Case I: Synchrotron spectrum with an exponential cutoff at high energy 

• Case II: Synchrotron spectrum with a broken powerlaw - aging population 

4.2 Modeling Implementation 

The radio through mid-IR data were modeled using a with a linear combination of the 
three models listed above (CLUMPY torus, starburst, and synchrotron emission). The 
modeled SED covers ~ 10 9 — 10 14 Hz (3/tm — 30 cm) in the rest frame of Cygnus A. The 
best fit combination (using the reduced x 2 ) was determined using Levenberg-Marquardt 
minimization. Some parameters are variables within the code and can be continuously 
varied while others must be changed by loading a new text file corresponding to that 
parameter set and model. 

An output from the CLUMPY code was loaded with a specific set of parameter val- 
ues. The same was done for the Siebenmorgen starburst models. An initial guess was 
made at the bolometric AGN flux (Fagn), the synchrotron cutoff frequency (z/ c ), and 
the optical depth of a dust slab attenuating the radio (r r ). The flux was calculated for 
these combination of parameters, and Levenberg-Marquardt least squares minimization 
was used to converge on the minimum x 2 value. The "best" values of Fagn, v c , and r r 
were recorded with the parameters from the loaded CLUMPY and starburst models, along 
with the reduced x 2 ■ Then, one of the parameters was varied in either the CLUMPY or 
starburst model (and the appropriate files loaded), and the process was repeated. The full 
list of parameters and their type (variable vs fixed) is given in table 14.21 If a parameter 
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is limited to discrete values, the range of values is given in the table. Otherwise "yes" is 
listed, indicating the variable is continuous. 

As there were a wide variety of parameters among the various models included, it was 
advantageous to constrain the parameters as much as possible prior to fitting. Some had 
constraints from previous studies or our data. The inclination of Cygnus A is known to be 
between 50° < i < 85°, constraining the choices of orientation for the CLUMPY model. 

The core radio fluxes specify the power-law index before the break as well as the zero 
point. A powerlaw was fit to the radio core fluxes and used to determine the spectral 
index (a = 0.18) and the flux scaling. 

The result of the combination of models was a grid of roughly 4 x 10 6 models. Due 
to computational constraints, it was not possible to fit each individual model combina- 
tion to the data. Accordingly the model grid was sampled in all parameters to yield a 
more manageable selection. The sampled grid consisted of approximately 842,000 model 
combinations. 
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Model 


Ta 

Parameter 


3le 4.2: Parameter List 
Type Continuous?/Range 


CLUMPY 


Fagn 

N 

Y 

q 

a 

TV 

i 

AGN PL visible? 


variable yes 

variable 1-25 

variable 5, 10, 30, 100, 200 

variable 0-5 

variable 10-85 (degrees, in 5° increments) 

variable 10, 20, 30, 40, 60, 80, 100, 150, 200, 300, 500 

variable 0-90 (degrees, in 10° increments) 

variable on /off 


Starburst 


r 

Ltot 
Job 
A v 
n 


variable 0.35, 1, 3, 9, 15 (kpc) 

variable 10 10 - 10 14 (L ) 

variable 0.4, 0.6, 0.9 

variable 2.2, 4.5, 7, 9, 18, 35, 70, 120 

variable 10 2 - 10 4 (cm -3 ) 


Synchrotron 


fo 
V Q 

1 v break 

ai 

Oil 
T r 


fixed NA 
fixed NA 
variable yes 
fixed NA 
fixed NA 
variable yes 
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Figure 4.3: \F\ vs A SED for the CLUMPY model, showing the effect of changing the 
various parameters on the output. The base parameters were: iV=20, Y=200, q=3, a=60, 
ry=40, accretion disk powerlaw not visible, and i=60. The following parameters are varied 
in plots, starting at the top: i, accretion disk powerlaw visible/not visible, a, q, Y, N, Ty . 
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Figure 4.4: Effect of model parameters on the output of the lSiebenmorgen &; Kriige 3 (|2007l ) 
models. The following parameters are varied with the constant values in parentheses. From 
top: A v (L tot = lO ia5 L , r=3 kpc, f OB = 0.9, n = 10 4 cm" 3 ); n (L tot = 1O 1O 5 L , r=3 
kpc, A v = 9, /on = 0.9); foB (Ltot = 1O 1L1 L , r=3 kpc, A v = 4.5,n = 10 4 cm" 3 ); L 



tot 



(r=3 kpc, A v = 17, f OB = 0.6, n = 10 3 cm" 3 ); r (L tot = 10 n i L Q , A v = 4.5, f OB = 0.9, 
n = 10 4 cm -3 ). The flux given is the spatially integrated flux at a distance of D[=50 Mpc. 
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Chapter 5 



Results 



The synthesis of physical models (each with multiple parameters) into a coherent picture 
of the physical nature of an object is unlikely to be as obvious as desired. A degree of 
degeneracy can be expected (and indeed is seen) when attempting to fit a ~15 parameter 
model. However, even with the ambiguity, the fits still shed light on the physical processes 
and conditions in the nucleus of Cygnus A. 

Figures 15.11 - 15.61 show the best 10 fits to the data for % = 80, 70, 50° with either an 
exponential cutoff (Case I) for the synchrotron spectrum or a break in the powerlaw (Case 
II). All fits shown have a x 2 /DOF < 2.3. The models accurately reproduce the radio, 
sub-mm, and IR emission seen in the nuclear region. 

Each collection of model fits managed to fit the data somewhat accurately. Table 15.11 
contains the average reduced x 2 value for the top 10 fits in each class. 

Table 5.1: Average reduced x 2 values for the top 10 fits in each set. 



i 


Case I 


Case II 


50 


1.969 


2.062 


70 


1.828 


1.317 


80 


2.266 


1.307 
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Figure 5.1: Cygnus A SED from radio through 4.5/im with the 10 best fit models for 
i = 80°, using a synchrotron spectrum with an exponential cutoff (Case I). Data from this 
work and references given in Ch[H 

5.1 Radio 

The source of the power-law emission in Cygnus A has been understood for some time 
as non-thermal synchrotron emission from a population of relativistic charged particles. 
This dominates from radio frequencies through several hundred GHz. At higher frequencies 
thermal emission begins to dominate the observed SED. 

5.2 Infrared 

The non-thermal synchrotron radiation obvious below a few hundred GHz contributes to 
the total FIR flux, but drops off quickly after ~10 THz, when the break in the particle 



53 



10 




0.00 1 0.D1 D.1 I 10 100 1000 

FiecpjenGV fTHi) 

Figure 5.2: Same as figure [5TT1 but for i = 70°. 

energy distribution is reached. However, even before the break, thermal emission from the 
torus and starburst dominates the synchrotron. 

The model components selected successfully reproduce the observed infrared emission 
(Figures 15.11 - I5T61) . The details of each component will be discussed below, focusing on 
the overall results and implications of the model components. § 15.2.11 discusses Case I fits 
(those with an exponential cutoff to the synchrotron population). § 15.2.21 shows results 
from Case II fits (those with a break in the powerlaw due to aging of the relativistic 
electron population). 

5.2.1 Synchrotron Properties and Contribution — Case I 

The synchrotron emission amplitude was fixed using the non-thermal emission from the 
radio core, assuming it to be a point source at all frequencies observed. As a result, there 
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Figure 5.3: Same as figure I5TTI but for i = 50°. 

were only two parameters fit: spectrum cutoff and extinction due to dust. None of the 
best fits had significant attenuation of the synchrotron due to dust. 

The AGN is behind a dust lane with significant extinction [Ay = 50±30, Ak = 5.6±3.3 



Djorgovski et al 



1991 



). Though some attenuation of the synchrotron flux is expected at 
higher frequencies, this alone is not sufficient to explain the cutoff. With extinction alone, 
the flux densities at 30 and 10/um would be similar for the spectral index given. This 
would exceed the measured infrared flux at 10/xm, so there must be a break in the electron 
population. 

In general, the cutoff frequency varies with the inclination angle of the torus, with 
some scatter for a given inclination. Edge on inclinations generally show a cutoff around 
10 THz (30^m). At this cutoff frequency, non-thermal emission contributes ~10% of the 
observed flux in the FIR and MIR. 
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Figure 5.4: Cygnus A SED from radio through 4.5/im with the 10 best fit models for 
i = 80°, using a synchrotron spectrum with a broken powerlaw resulting from aging of the 
relativistic electron population without particle injection (Case II). Data from this work 
and references given in Ch|H 

Less extreme inclinations (i ~50-70) show slightly higher cutoff frequencies, on the 
order of 18-20 THz (15-16 fim). At these inclinations, the synchrotron again contributes 
~10% of the observed flux in the FIR and MIR. However, at frequencies higher than the 
cutoff frequency, the synchrotron can contribute over 50% of the observed flux (for exam- 
ple, around 30 THz, 10/xm). With such a significant synchrotron component, variability 
around 10/x m might be present. 

Assuming each electron emits only at its critical frequency, the gamma factor for an 
electron is given by eq 15.11 (where B is the magnetic field strength, v is the frequency). 

/ 4-Km e cu 

7 = V^~ (5 - 1} 
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Figure 5.5: Same as figure [5T41 but for i = 70°. 

Plugging in numerical values, this becomes eq 15.21 (for v in THz and B in /i Gauss. 
The energy of the electron can then be computed using E = jmc 2 . 



7 = 4.884 x 10 s 



v \\ ( B \~t 



THz 



(5.2) 



Table 5.2: Synchrotron Cutoff Frequencies 



i 


Median v c 


Median A c 




^/cutoff 


Lsync 


deg 


THz 


[xm 




B=100^G 


v ~ 10 9 - W u Hz; erg s" 1 


80 


10.62 


28.05 


5.46 


1.59 x 10 a 


2 x 10 44 


70 


19.84 


15.02 


0.14 


2.17 x 10 5 


4 x 10 44 


50 


19.19 


15.52 


0.42 


2.14 x 10 5 


4 x 10 44 



There is clearly signifi c ant ex tinction in the optical due to the dust lane across the 



nucleus. iDiorgovski et al 



(|199ll ) derive Ak = 5.6 db 3.3 for the nucleus. One possible 
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Figure 5.6: Same as figure [53] but for i = 50°. 

explanation for the break in the synchrotron flux is simply absorption due to the dust lane. 
Figure O (LEFT) shows the opacity scaled to A K = 5.6. Figure O (RIGHT) plots a 
synchrotron spectrum with a = 0.18 behind a dust lane with Ak = 5.6. Comparison of the 
attenuated powerlaw with the emission around 50 THz indicates absorption alone cannot 
account for the observed spectrum of Cygnus A. The powerlaw must break somewhere in 
the mid- infrared. 

The best fit synchrotron parameters for a powerlaw with an exponential cutoff are 



listed in table 15.21 7 is calculate d for a 100 nG fie 



are ~ 130 fiG in knots in the jet (ICarilli &: Barthel 



d. M easured values of magnetic field 



1996). These calculations assume the 



"minimum energy" condition where the energy in relaivistic particles is roughly equivalent 
to the energy in magnetic fields. 100 fiG was adopted as a round number on the order 
of the magnetic field strength in the kpc-scale jet. Note that magnetic field strengths 



58 



measured in the VLBI radio core are higher, 0.16 G. The lower value measured using VLA 
observations was chosen as the VLA "core" includes emission from both the VLBI core 
and the jet. 
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Figure 5.7: Lef t : Th e braine Lei opacity curve scaled to Ak given by 

Djorgovski et alJ (jl99ll ) (dashed line). Also plotted is a F u oc v~ 0A8 e~ T (solid line), 
representing the measured synchrotron spectrum of the core attenuated by a dust slab. 
Right: The same powerlaw, overplotted with the SED of Cygnus A. An absorbing slab 
consistent with that is seen in Cygnus A is not sufficient to explain the observed properties 
of the synchrotron emission. 



Energy loss in synchrotron is proportional to both the square of the magnetic field 
strength and the square of the particle energy. For a given particle energy, a synchrotron 
cooling time can be defined as in eq[ 



T. 



sync 



( B 
635 - 

v< 



-)" 

erg J 



(5.3) 



seconds \ Gauss 

For the magnetic field assumed above (100^G) and a particle energy at the cutoff 



(given by 7 above), typical synchrotron cooling times are roughly 12,300 years. 
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Inverse Compton Scattering 

In addition to producing synchrotron radiation, a population of relativistic electrons will 
scatter photons to higher energies through inverse compton scattering with power given 
in eq 15,41 ot is the Thomson cross section, f3 = v/c and U ra d is the energy density in the 
radiation field. 

Pic = ^ T cp 2 7 2 U rad (5.4) 

For comparison, synchrotron losses for relativistic electrons follow a similar relation, 
with the radiation energy density replaced by the magnetic energy density (eq 15.51 where 
Ub is the energy density in the magnetic fields). 

Psyn = ^TC(3 2 7 2 U B (5.5) 

Thus, the ratio of the energy lost through synchrotron to the energy lost to IC is 
simply given by the ratio of the energy densities. The Cosmic Microwave Background 
(CMB) sets a minimum value for the radiation energy density. Using the magnetic field 
strength assumed above (100// G), and Tqmb = 2.73K for the CMB, synchrotron losses 
are dominant by a factor of ~ 1000. A larger magnetic field would further enhance this, 
so IC losses from scattering of CMB photons have a negligible effect on the population of 
relativistic particles. 

The resulting frequency of IC scattered CMB photons is proportional to j 2 . A Tqmb = 
2.73K blackbody has an emission peak at vcmb ~ 160 GHz. CMB photons at the peak 
which scatter off relativistic electrons emitting at the synchrotron cutoff frequency would 
be scattered to frequencies of vjc ~ 8.5 x 10 21 Hz (from eq 15.61) . This corresponds to 
photon energies of ~ 35 MeV. IC scattering off a powerlaw distribution of relativistic 
electrons (as in Cygnus A) results in IC emission with the same spectral index (a = 0.18). 
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VIC = ^l 2 Vrad (5.6) 

In addition to scattering of CMB photons, the population of relativistic electrons can 
interact with the synchrotron radiation via IC in what is known as "Synchrotron Self- 
Compton" (SSC). Synchrotron radiation at the cutoff frequency would be scattered to 
frequencies of v ~ 6 x 10 23 Hz [E ~ 2.4 GeV). 

5.2.2 Synchrotron Properties and Contribution — Case II 

An alternate mechanism for limiting the influence of synchrotron emission at shorter wave- 
lengths is for the spectrum to break at some frequency (see § I4.1.3P . An identical process 
to that above was used to fit a broken powerlaw synchrotron spectrum (plus torus and 
starburst models) to the Cygnus A data. 

Fits used a pre-break spectral index of a\ = 0.18 and a post-break spectral index 
of «2 = 1-24, consi stent with aging of the relativistic population (without injection of 



additional particles; 



Kardashev 



19621 ) . As the flux density decreases in a slower fashion 
when compared to an exponential cutoff, the powerlaw must break at lower frequencies to 
ensure the observed midIR flux is not exceeded. Table 15.31 is analogous to Table 15.21 but 
for the broken powerlaw. 

Table 5.3: Synchrotron Break Frequencies 



i 


Median u break 


^break 


break 




Lsync 


deg 


THz 


/j,m 




B=100/xG v r 


- 10 9 - W u Hz; erg s" 1 


80 


7.44 


40.04 


0.42 


1.33 x 10 5 


5 x 10 44 


70 


6.91 


43.11 


0.33 


1.28 x 10 5 


5 x 10 44 


50 


3.16 


94.28 


0.55 


0.868 x 10 5 


3 x 10 44 
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Inverse Compton Scattering 

Similarly to Case I, IC from a relativistic electron emitting at ~ 5 THz will be scattered 
to vie ~ 8.0 x 10 22 Hz (E ~ 330 MeV). Potentially higher energy photons could result 
from scattering off electrons in the post-break distribution. 

5.2.3 AGN/Torus Properties and Contribution 

An estimate for the inclination of the AGN was given based on previous work in the 
radio regime. This constraint was used in fitting the SED. However, acceptable fits were 
obtained for all inclinations within the range, so no further constraint can be placed on 
i using this analysis. For the each of the inclinations, the properties of the best fit torus 
model will be discussed below. Summary tables are provided for both Case I and Case II 
fits following the discussion of the individual models (Tables 15.41 & I5.5f) . 

i=80° 

Figure [5TT1 shows the Cygnus A SED with the 10 best fit Case I models at i = 80. These 
best fits are individually plotted in GTS All models had 1.802 < x 2 /DOF < 1.865. In 
general, the best fit torus model did not have contributions from the AGN powerlaw. 

For model combinations with the AGN directly visible, the bolometric AGN flux is 
~ 680 x 10 -10 erg s _1 cnr 2 . At a luminosity distance of 237 Mpc (7.31 x 10 26 cm) the 
luminosity of the AGN is L agn = 4.6 x 10 47 erg s _1 . Without the AGN directly visible, 
the AGN flux is ~ 14 x 10 -10 erg s _1 cm -2 , giving L agn = 9.4 x 10 45 erg s _1 . The inner 
radius of the torus is then 8.6 pc for the fits with the AGN powerlaw visible and 1.22 pc 
for fits where the powerlaw isn't visible. 

The CLUMPY torus model which best fits the data at i = 80 includes a visible AGN 
powerlaw. Due to the clumpy nature of the torus, there is a finite probability of directly 
viewing the AGN. These fits involve a relatively thin torus (a = 15) with a large number 
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of clumps along the line of sight (N = 20), each with moderately high optical depth 
(t v = 40). The torus is also quite broad in the radial direction (Y = 200) and the density 
of clumps does not drop very rapidly. 

These models with the AGN visible require a much more luminous AGN (by at least 
one order of magnitude, sometimes two). The bulk of the energy output (> 90%) in these 
fits is in the optical and UV, which is not constrained here. Given such optical AGN 
continuum is not observed, these fits do not seem physically plausable for Cygnus A. 

Mid-infrared observations of NGC 1068 (also a type II AGN) indicate the torus is 



contained within a 15 pc region ([Mason et al 



20061 ). Fits at i = 80 with the AGN pow- 



erlaw visible require an outer torus radius for Cygnus A of ~ 1.7 kpc! This is unlikely 
to be a realistic result, as a 1.7 kpc torus would have easily been observed (at the dis- 
tance of Cygnus A, 1" corresponds to ~ 1 kpc, so a structure of this size would be well 
within range of HST and some ground-based observatories). Additionally, these fits re- 
quire a bolo metric AGN luminosity at least a factor of 10 larger than implied by x-ray 



observations ( Tadhunter et al 



2003) , and a factor of over 100 



luminosity of the AGN implied by 



Whysong &: Antonuccil (|2004l ) . Finally, the torus model 



arger than the bolometric 



component for t 



lese fits significantly overshoots the nuclear K-band spectra presented by 



Tadhunter et al 



(|2003h 



For the model in which the AGN is not visible, the torus model was markedly different 
in physical configuration, being much smaller (Y = 30; a = 60), but having otherwise 
similar properties. 

The NIR provides an opportunity to break the degeneracy between models with the 
AGN powerlaw visible and those without the AGN powerlaw. The AGN powerlaw results 
in an upturn in the SED around 100 THz (3/im), while other models have a flux den- 
sity which continues to fall. Preliminary fitting using 2MASS points indicates the AGN 
powerlaw exceeds the observed flux in the NIR. 
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Figure 15.41 shows the top 10 Case II fits at i = 80. These models had 1.273 < 
X 2 /DOF < 1.336. The torus model fits are very similar to the Case I fits where the 
AGN powerlaw is not visible, with an AGN that is ~ 15% more luminous. Individual 
plots of the fits are in figure 15.91 

i=70° 

At the slightly shallower inclination of i = 70 a much less luminous AGN is required to fit 
the data (Case I fits shown in figures IST21 and [5,10p . Typical bolometric AGN fluxes were 
~ 1.15 x 10 -10 erg s -1 cm -2 , corresponding to a luminosity of ~ 7.7 x 10 44 erg s . x 2 
values were ~ 1.8. The corresponding inner radius of the torus is ~ 0.35 pc. Reduced x 2 
values are between 1.780 and 1.858. 

As with i = 80, torus models with and without the AGN powerlaw provide acceptable 
fits to the available data. However, in contrast, they both predict similar bolometric 
fluxes for i = 70. The best fit torus models have a relatively small opening angle, with 
a ~ 60 — 65°. The other fit parameters are Y = 30, ./V = 22, q = 2 and t v = 60. 

Case II fits (figures 15.51 & 15. 12j) select torus models which are effectively identical 
to those in Case I fits in the physical characteristics. The only difference occurs in the 
luminosities, which tend to be ~ 50% higher in the Case II fits. Reduced x 2 values are 
between 1.290 and 1.345. 

i=50° 

Case I fits using an inclination angle of i = 50 imply a bolometric AGN flux similar to 
that of i = 70: ~ (0.85 — 1.6) x 10~ 10 erg s _1 cm~ 2 , corresponding to a luminosity of 
(5.8 — 11) x 10 44 erg s _1 . The favored torus model is fairly consistent among the best fits. 
At this luminosity the inner radius of the torus is ~ 0.3 — 0.4 pc. The reduced x 2 values 
are between 1.888 and 2.006. 
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The best fit torus parameters for i = 50 are similar to those for i = 70, with slight 
fluctuation in the powerlaw distribution of clumps (q between 1 and 2), and a slightly 
wider range in the angular distribution of clumps (a between 55 and 65°). As with the 
other inclinations, Case II fits result in similar torus properties as Case I fits for i = 50 
(Figures 15.61 &: [5.14p . In this case the best fit luminosities are similar. Reduced x 2 values 
are between 1.997 and 2.122 for the top ten fits. 

The best fit torus parameters for the inclinations are summarized for Case I in Table 
15.41 and Case II in Table 15.51 

Table 5.4: Summary of Torus Model Results: Case I - Exponential Cutoff 



i 


AGN? 


Y 


a 


N 


q 


r v 


Lagn (erg S" 1 ) 


V = LAGN/L E dd 


80° 


yes 


200 


15° 


20 


1 


40 


4.6 x 10 4Y 


0.56 


80° 


no 


30 


60° 


20 


1 


40 


9.4 x 10 45 


1.1 x 10" 2 


70° 


yes & no 


30 


60-65° 


~20 


2 


60 


7.7 x 10 44 


9.3 x 10~ 4 


50° 


yes & no 


30 


55-65° 


~20 


2 


60 


(5.8-11) xlO 44 


7.0 - 13. x 10~ 4 



Table 5.5: Summary of Torus Model Results: Case II - Broken Powerlaw 



i 


AGN? 


Y 


a 


JV 


q 


T v 


Lagn (erg S _1 ) 


V = LAGN/LEdd 


80° 


yes & no 


30 


60° 


20 


l 


40 


1.1 x 10 45 


1.3 x 10~ 3 


70° 


yes & no 


30 


60° 


20 


l 


60 


1.1 x 10 45 


1.3 x 10~ 3 


50° 


yes & no 


30 


60° 


20 


l 


60 


7.3 x 10 44 


8.8 x 10~ 4 



Physical Interpretation 



Clearly the best fit models have settled on a torus with a small opening angle and a 
shallow falloff in clump density with radius. The optical depth of each clump is also 
fairly high, but not excessively so. Additionally, with the exception of i = 80 + visible 
AGN, the luminosities are similar to the bolometric AGN luminos ity inferred from X-ray 



and mid-IR observations (gi ving = 0.1 — 2 x 10 ergs s 1 ) (jTadhunter et al 



Whysong k, Antonucci 



200J) 



2003 
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Young et al.l (|2002l ) found Nh = 2.0 x 10 23 cm~ 3 towards the nucleus from X-ray 
observations. For a standard gas-to-dust ratio, this column density corresponds to Ay ~ 
100 mag. Torus model fits find a somewhat lower Ay ~ 65. 

There is little difference in the Case I and Case II fits with respect to torus models. 
Torus emission is dominant in the midIR, where the synchrotron arguably has the smallest 
effect in either case. 



5.2.4 Starburst Properties and Contribution 

The best fit starburst component varies significantly between the inclinations, however 
remains relatively consistent for the va rious fits a t a give n inclination. A star formation 



rate is calcul ated using Lftr following 



estimates see iKennicuttl (|1998al ) 



Kennicuttl (|l998bh (eqEZD- For a review of SFR 



SFR (M yr' 1 ) = 4.5 x 10 



-44 LpiR 



1 erg s 1 



1.72 x 10" 



-10 



Lfir 



(5.7) 



At nearly edge on inclinations i = 80, a starburst of L s b ~ 10 L & with A v ~ 18 is 
the best fit (with scatter in A v from 4.5 through 72). The size is typically 3kpc (ranging 
from 0.35-3 kpc), similar to the size of the circumnuclear ring known to exist in Cygnus 
A (4kpc). The dust density in hot spots is typically n ~ 5 x 10 3 cm 3 , with OB stars 
contributing roughly 60% of the luminosity. Using eq 15.71 the SFR rate is ~27 M yr _1 . 

With i = 70, the best fit starburst model changes slightly, favoring a smaller burst 
(r = 0.35 kpc), and a slightly lower Job = 40%. The total luminosity decreases somewhat 
when compared to fits for an torus orientation of i = 80: L s b ~ 10 11 L Q , while the 
extinction increases to A v = 35 (scatter from 18-72). The dust density in the hotspots 
ranges from 1 — 10. x 10 3 cm -3 , and the SFR is 17 M yr _1 . 

The fits with a torus inclination of % = 50 are very similar to those of i = 70 with a 
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smaller range in the density in the hotspots of the OB stars (1 — 2.5 x 10 3 cm -3 . The 
luminosity is also slightly higher: L s b ~ 10 L©, corresponding to a SFR of 22 M yr _1 . 

In contrast with the torus results, the best fit starburst properties do change somewhat 
between Case I and Case II fits. 

Table [5U1 summarizes the results of the starburst fits for Case I. Table [57T1 corresponds 
to Case II fits. 

Table 5.6: Summary of Starburst Model Results for Case I 



i 


r (kpc) 


foB 


Ltot (£©) 


A v 


n (cm 3 ) 


SFR (Meyr- 1 ) 


80 


3 


60% 


1.5 x 10 11 


18 


5 x 10 3 


27 


70 


0.35 


40% 


0.8 x 10 11 


35 


2.5 x 10 3 


17 


50 


0.35 


40% 


0.8 x 10 11 


35 


1 x 10 3 


22 



Table 5.7: Summary of Starburst Model Results for Case II 



i 


r (kpc) 


foB 


Ltot (Lq) 


A v 


n (cm 3 ) 


SFR {M & yr- 1 ) 


80 


0.35 


40-90% 


1 x 10 11 


72 


1 x 10 2 


17 


70 


0.35 


40-90% 


(0.8 - 1.7) x 10 11 


72 


~ 5 x 10 3 


14-30 


50 


3 


40% 


2 x 10 11 


4.5 


5 x 10 3 


34 



Physical Interpretation 

The best-fit starburst models at all inclinations are very similar in most areas. Roughly 
half the luminosity is from OB stars enshrouded in hot spots with dust densities of a 
fewxlO 3 cm -3 and having a total luminosity of ~ 10 11 L© (3 x 10 44 erg s _1 ). Star 
formation rates are roughly the same in all fits, ranging between 20 and 30 M© yr~ l for 
Case I and Case II fits. 

The integrated far infrared luminosities of the starburst component gives significant 
star formation rates, indicating Cygnus A would still be bright in the infrared, even in the 
absence of an AGN. The luminosities are comparable to those found in luminous infrared 
galaxies (LIRGs). 
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5.2.5 Far IR Analysis 

Around 1 THz (~ 300//m) the three model components are of similar flux density, with the 
synchrotron dying out and the thermal contributions from star formation and a clumpy 
torus strengthening. Unfortunately, data at this location is quite limited, and of relatively 
poor spatial resolution, leaving little hope of spatially disentangling the emission from 
the starburst and emission from the AGN at this frequen cy until the availability of the 



Atacama Large Millimeter Array (ALMA; eg 



Carilli 



2005J) 



The thermal emission seen in the far-IR dominates at frequencies higher than 1 THz 
(~ 300/iin). The reprocessed emission from the AGN exceeds the starburst emission just 
shy of 10 THz (~ 30/um). At this point, the flux density of the observed SED is decreasing, 
and the peak flux density of the torus contribution is lower than the peak flux density of 
the starburst. 

5.2.6 Mid IR 

The MIR portion of the spectrum contains many (dust) emission features which are ex- 
pressed in the models. Additionally, there are numerous forbidden lines from highly ionized 
species as well as H2 molecular lines in the IRS spectrum. These are likely from the inner 
NLR, not the AGN/torus region (BLR), and discussion of them is beyond the scope of 
this study. 

An obvious feature of the MIR spectra is the broad absorption around lOfim. This 
silicate absorption is well matched by the combination of torus and starburst models (see 
inserts in figs 15.81 - 15.14ft . 

As noted earlier, the continuum flux in this regime requires the synchrotron powerlaw 
to break at an unknown frequency. Based on the above modeling results, the break 
likely occurs in the mid- infrared (~ 15/um). In cont rast to Cygnus A, some AGN exhibit 



optical synchrotron jets (eg. 3C133, 3C273, M87; 
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Floyd et al 



2006; 



Jester et al. 



2007 



Perlman et al.ll2007l ). indicating the distribution of the relativistic electrons might extend 
to higher energies in these sources. 

Approaching the NIR, all model components continue to rapidly decrease in flux den- 
sity. In fits where the AGN powerlaw is visible, this component begins to dominate into 
the NIR. However, the stellar contribution to the NIR have not been taken into account. 



5.3 Overall Contributions to the Luminosity in Cygnus A 

The three physical processes included in modeling the radio through mid-IR SED of 
Cygnus A all contribute to the luminosity. Of the three, torus reprocessed AGN radi- 
ation dominates the energy output of the nuclear region. Table 15.81 lists the ratios of each 
model component. 

Table 5.8: Ratio of Luminosity for Model Components 





i 


Lagn 


Lsb/Lagn 


L S ync/ Lagn 




80 


4.6 x 10 4Y 


1.2 x 10~ 3 


8.7 x 10" 4 


Case I 


80 


9.4 x 10 45 


0.061 


0.043 


70 


7.7 x 10 44 


0.39 


0.52 




50 


~ 8 x 10 44 


0.37 


0.25 




80 


1.1 x 10 45 


0.38 


0.27 


Case II 


70 


1.1 x 10 45 


0.43 


0.45 




50 


7.3 x 10 44 


1.0 


0.68 



In most cases, the re-radiated emission from the torus dominates the luminosity by a 
factor of a few. In addition, the energy emitted in the core through synchrotron is roughly 
the order of that in the starburst, and factors of a few less than the emission from the 
torus. 
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0.00 1 0.01 D.1 I 10 10D 0.001 D01 0.1 I 10 100 



FreqdSncy (THz) 

Figure 5.8: Individual plots of the best fit models in figure [5TT1 Legend is the same 
figure I5.ll 




FraquQncy (T Hz) 
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0.00 1 0.01 D.1 I 10 100 0.DQ1 DQ1 0.1 I 10 100 



Frequency (THz) 

Figure 5.9: Individual plots of the best fit models in figure [53 Legend is the same 
figure I5.4I 




FraquQncy (T Hz) 
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0.00 1 0.01 D.1 I 10 100 0.K31 001 0.1 I 10 100 



Frequency (THz) 

Figure 5.10: Individual plots of the best fit models in figure [531 Legend is the same as 
figure 15.21 




djgdi a. 01 d.1 1 ia iaD a.wi o.ai a.1 i 10 100 



FnaquQncy (THz) 

Figure 5.11: Individual plots of the best fit models in figure I5T21 Legend is the same as 
figure 15.21 
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0.00 1 0.01 D.1 I 10 10D 0.D01 DCM 0.1 I 10 100 



FreqillVicy (THz) 



Figure 5.12: Individual plots of the best fit models in figure [531 Legend is the same as 
figure 15.51 




FraquQncy (T Hz) 
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a. oi 



D.DQ1 - 




□.□1 - 



D.H31 - 



0.001 



o.oo i a. ai oi 



10 100 0.D01 001 0.1 I 10 100 



Freqil^ncy (THz) 



Figure 5.13: Individual plots of the best fit models in figure [5731 Legend is the same as 
figure [ 




FraquQncy (T Hz) 
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0.00 1 0.01 D.1 I 10 100 0.K31 001 0.1 I 10 100 



FreqilrSncy (THz) 

Figure 5.14: Individual plots of the best fit models in figure [5U1 Legend is the same as 
figure 15.61 




FraquQncy (THz) 
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Chapter 6 



Conclusions 



Using a combination of new mid-infrared spectra from the Spitzer Space Telescope and 
radio data from the literature, the nuclear emission of Cygnus A has been modeled as a 
combination of powerlaw emission from a synchrotron jet, reprocessed AGN emission from 
a dusty torus, and reprocessed emission from a circumnuclear starburst. This is one of 
the first detailed fits to a broad SED, covering ~ 5 dex in frequency. The addition of the 
mid-IR data is critical in constraining the relative contributions of each component to the 
bolometric luminosity of Cygnus A. It places limits on the amount of warm dust emission 
present in the inner regions, constraining the contribution from the torus component. 

The data is well fit by a combination of these three models. The fits suggest the nuclear 
emission at these frequencies is primarily due to AGN activity (synchrotron radiation & 
torus reprocessed AGN emission) and star formation. 

The modeled jet shows evidence for a cutoff in the mid-infrared. This is in contrast 
with some radio galaxies and quasars which showing jets emitting from the radio through 
X-rays. Statistically acceptable fits were found for both an exponential cutoff in the 
population of relativistic electrons (Case I) as well as synchrotron losses (Case II). The 
cutoff/break frequencies tended to be within the 3 — 30 THz range. It is possible the 
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synchrotron emission is a significant component of the mid-IR flux, which could be tested 
using polarimetry. 

The nuclear AGN emission was modeled as a clumpy torus, reprocessing the emission 
from and around the accretion disk. The fits suggest a bolometric AGN luminosity of 
~ 10 45 erg s _1 . The torus is predicted to have an outer radius of ~ 7 pc (7 mas at the 
distance of Cygnus A). 

The star formation rate determined for Cygnus A is on the order of 10 — 30 M yr _1 , 
well above the rate for quiescent galaxies. The IR luminosity of this component is large 
enough that this galaxy would be classified as a LIRG, even in the absence of AGN activity. 

The reprocessed torus radiation dominates the bolometric luminosity of Cygnus A, 
typically a factor of ~ 3 greater than the luminosity associated with star formation. While 
the dominant emission at radio wavelengths, the synchrotron emission is a small fraction 
of the emission in the far-IR. 

It was noted that Cygnus A is in the midst of a merger. The existence of the merging 
system, star formation, and AGN activity is consistent with suggestions that mergers 
funnel gas towards the nucleus where it forms stars and feeds a SMBH at the center. 

6.1 Further Work 

There remains some ambiguity in the fits primarily related to the visibility of the AGN 
powerlaw. Inclusion of data in the near infrared should resolve this issue. Preliminary use 
of 2MASS data indicates the AGN powerlaw is not visible. However, this remains to be 
properly verified by taking into account stellar contribution in the near-infrared. 

The contribution of synchrotron self-compton (SSC) to the observed X-ray flux remains 
unquantified. Synchrotron photons scattering off the relativistic particles that emitted 
them can potentially be a significant contribution to the X-ray luminosity in Cygnus A. 
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The Herschel Space Observatory, set to launch in May 2009, will provide spectral 
and imaging coverage in the far infrared and sub-mm. Low- and medium- resolution spec- 
troscopy (i? ~ 40—5000, depending on instrument and wavelength) will cov er a wavelength 
range of 55 — 672/im, pr obing a relatively unexplored wavelength regime (jPoglitsch et al. 



2006 



Griffin et al 



2006). This will likely put additional constraints on the physical pro- 
cesses contributing in the far infrared. Also, atomic and molecular lines at the longer 
wavelengths will provide additional information on the physical conditions present in the 
nuclear regions. 

The sub-mm is a wavelength regime in which there is a comparable dearth of informa- 
tion. Historically plagued with poor angular resolution, future instruments, particularly 
ALMA will enable deeper, higher-resolution studies of the nuclear regions. It will be able 
to resolve and image the nuclear region, potentially spatially disentangling the circumnu- 
clear starburst and AGN, removing some ambiguity in the modeling. However, Cygnus A 
is at a high declination and might prove difficult to observe with ALMA. 

The launch of the James Webb Space Telescope (JWST) will facilitate deeper and more 
sensitive observations in the near- and mid-infrared. The Mid-Infrared Instrument (MIRI) 



will cover a wavelength range similar to tha t of IRS whi 



(both spatial and spectral) and sensitivity (jWright et al 



e offering increased resolution 



2004). This would enable a 



study of the nuclear region which resolves physical separation between the starburst ring 
and the AGN, testing the conclusions from the model fitting shown here. 
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